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Abstract 

If neutrino masses have a radiative origin, their smahness can be naturahy 
understood even when lepton number violation occurs near the weak scale. 
We analyze a specific model of this type wherein the neutrino masses arise 
as two-loop radiative corrections. We show that the model admits the near 
bimaximal mixing pattern suggested by the current neutrino oscillation data. 
Unlike the conventional seesaw models, these two-loop models can be directly 
tested in lepton flavor violating decays r ^ 3/U and ^ ^ e + 7 as well as at 
colliders by the direct observation of charged scalars needed for the mass 
generation. It is shown that consistency with the neutrino oscillation data 
requires that the leptonic rare decays should be within reach of forthcoming 
experiments and that the charged scalars are likely to be within reach of the 
LHC. 
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I. INTRODUCTION 



Recent neutrino oscillation experiments, notably Super-Kamiokande and SNO, have pro- 
vided a wealth of information on neutrino masses and mixing pattern. First and foremost, 
the solar [1] and the atmospheric neutrino data [2] are incompatible with the hypothesis of 
zero neutrino mass. The inferred values of neutrino masses arc in the sub-eV range. The 
mixing angle 612 relevant for solar neutrino oscillations appears to be large (tan^ 612 ~ 0.4) 
[3], while the one for atmospheric neutrino oscillations ^23 ~ 7r/4 is near maximal. Re- 
actor neutrino experiments [4] constrain the third mixing angle ^13 to be relatively small, 
^13 < 0.16. It will be of great interest to seek a theoretical understanding of these observa- 
tions. 

The seesaw mechanism [5] has been the most popular explanation of the small neutrino 
masses. It is elegant and simple, and relies only on dimensional analysis for the required 
new physics. Current neutrino data points to a seesaw scale of Mji ~ 10^° — 10^^ GeV, 
where lepton number violation occurs through the Majorana masses of the right-handed 
neutrinos. With such a high scale, effects of lepton flavor violation in processes other than 
neutrino oscillation itself becomes extremely small. For example, the branching ratio for the 
decay — > e -|- 7 is of order 10~^° within the seesaw extension of the Standard Model. 

An alternative to the seesaw mechanism which also explains the smallness of neutrino 
masses naturally is the radiative mass generation mechanism [6-14]. In this approach, 
neutrino masses are zero at the tree-level, and are induced only as finite radiative corrections. 
These radiative corrections are typically proportional to the square of the charged lepton (or 
quark) masses divided by the scale of new physics A. The neutrino masses can then be in 
the sub-eV range even for A ~ TeV. In this case, lepton flavor violation in processes other 
than neutrino oscillations may become experimentally accessible. 

In this paper we propose to analyze in detail the experimental consequences of a spe- 
cific model wherein the neutrino masses arise as two-loop radiative corrections [8,9]. The 
magnitude of these induced masses are of order uiy ~ [(/^/i)/(167r^)^](m^/A), where /, h 
arc dimensionless Yukawa couplings and A is the scale of new physics. For / ~ /i ~ 0.1 
and A ~ 1 TeV, the neutrino mass will be of order 0.1 eV, which is compatible with the 
oscillation data. Within this model, it is required that f,h> 0.1 and A < 1 TeV, or else 
the induced neutrino mass will be too small to be relevant experimentally. These limits, 
along with the requirement of large solar and atmospheric neutrino oscillation angles, imply 
that lepton flavor violation in processes such as r — > 3/i and — > e-1-7 cannot be suppressed 
arbitrarily. We shall see that these decays are within reach of the next round of rare decay 
experiments. Furthermore, the upper limit on A implies that it is very likely that the new 
scalars predicted by the model for neutrino mass generation will be accessible for direct 
observation at the LHC and perhaps even at Run 11 of the Tevatron. 

A variety of radiative neutrino mass models exist in the literature. The Zee model of 
neutrino mass [6] is a popular example. In this model, neutrino masses arise as one-loop 
radiative corrections. This model appears to be incompatible with the large mixing angle 
MSW solution for the solar neutrino problem as it predicts the relevant mixing angle to to 
very close to 7r/4 [15]. In any case, one would expect rare leptonic processes to be more 
suppressed in the Zee model compared to the two-loop neutrino mass model [8,9]. This 
is because the neutrino masses in the Zee model, being of one-loop origin, are of order 
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"^i^ ~ [/fl'/(167r^)](m^/A), and in order for it to explain the atmospheric data, one must 
choose either / < 10^'^ if A < TeV, or A > 10^ GeV if / ~ 1. This should be contrasted 
with the two-loop neutrino mass model where the Yukawa couplings are needed to be order 
0.1 and A of order TeV. 

Another class of neutrino mass models that has been widely studied is supersymmetric 
models with explicit i?-parity violation [10]. In this case, it turns out that only one neutrino 
acquires a tree-level mass, the other two obtain masses as one-loop radiative corrections. 
The phenomenology of such models, especially when i?-parity violation is soft, arising only 
through bilinear terms, has been thoroughly investigated [13]. 

In Ref. [14] a general effective operator approach for small Majorana neutrino masses has 
been presented. It follows from that analysis that if Icpton number violation resides only in 
the leptonic sector and does not involve the quarks, then the only interesting model for the 
current neutrino oscillation data that also predicts observable lepton flavor violation signals 
is the two-loop model of Ref. [8,9]. This is the main reason for revisiting the two-loop model 
in the present paper. 

The outline of the paper is as follows. In Section II we will have a quick review of the 
two-loop neutrino mass model. In Section III we analyze quantitatively the constraints on 
the model arising from neutrino oscillation data. Here we also show how the model can 
accommodate the near bi-maximal mixing pattern preferred by the current data. In Section 
IV, we analyze the constraints arising from rare lepton decays. Here we present bounds on 
the masses and couplings of the charged scalars present in the model. Section V is devoted 
to the collider implications of the model. We provide our conclusions in Section VI. 

II. DESCRIPTION OF THE MODEL 

We start by reviewing the essential features of the two-loop neutrino mass model [9,8]. 
The gauge group of the model is the same as the Standard Model (SM). In addition to the 
SM particles the model introduces two charged scalars h'^ and k~^~^ which are both singlets 
of SU{3)c' and SU{2)l. These scalars have Yukawa couplings to the leptons, which can be 
parametrized in terms of two complex matrices / and h: 

C^Y = /ab ifaiCi^ij) ^ijh-' + K, {llRClm) + h.c. (1) 

Here V'l stands for the left-handed lepton doublet, and Ir for the right-handed lepton singlet. 
C is the charge conjugation matrix. i,j are SU{2)l indices, while a, h are generation indices. 
The matrix / is antisymmetric {fab — —fba) due to Fermi statistics and antisymmetry in 
the SU{2)l indices, while h' is a symmetric matrix {hl^^ — h'^g)- The interaction Lagrangian 
in terms of the component fields will then be 

Cy = 2 [/e^(i7c/Xi - V-^tL) + feMTL - U^Cl) + UA'^f.TL - U^PLl) ] (2) 
+ [ heee^CR + h^^fl^HR + hrrT^TR + he/.e^fXR + here^TR + hi^rf^^TR ] k'^^ + h.C. 

where we have defined haa = h'^a^ hah = 2/;,^^ for a 7^ 6. 

The Yukawa interactions of Eq. (1) conserves lepton number (L), as can be seen by 
assigning two units of L to the /i+ and A;"*"+ fields. The scalar piece of the Lagrangian 
contains a term 



3 



FIG. 1. Diagram leading to finite neutrino mass at two loop. 

Lh-}^ — —fjLh'^h'^k + h.c. (3) 

which would then violate L. In fact, one sees that the combination of Eqs. (1) and (3) 
explicitly breaks lepton number. This would lead to the generation of Majorana neutrino 
masses at the two-loop level. The relevant diagram is shown in Fig. 1. The induced neutrino 
mass can be calculated to be 

{M^)ab = ^l^facrricKai^dfdblcd (4) 

where rric, rud are the charged lepton masses and led is the two loop integral function given 
by 



d^k r d^q 1 1 1 1 . 

^-d- I (2^y (2^ {k-i - ml) {k^ - ml) {q^ - mj) {q^ - ml) {k - qf - ml ' ^ ' 



In order to evaluate the above integral, we may neglect the lepton massess in the denom- 
inator, since these massess are much smaller than the charged scalar masses m^ and m^. 
Then 

(167r2)2 ml \mlJ 
The dimensionless quantity / defined as 

/(.) = -f\4'-' iy^i^y i„/Ji^ (7) 



JO 



x+{r — l)y + y^ x + ry 



is plotted in Fig. 2 for r = ml/ml < 100. The asymptotic behavior for small and 
mk mfi is given by: 

I{r) ^ { (8) 

7rV3 for r ^ 

Before we estimate the induced neutrino masses and mixings, we wish to make some 
remarks on the expected magnitude of the relevant coupling parameters fab, hob and First, 
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FIG. 2. The dimensionless function I{r) appearing in the evaluation of the two-loop mass 
integral. 

in order for perturbative expansion to be sensible, the Yukawa couplings fab, hab have to be 
of order unity or smaller. For concreteness we shall take fab < 1, hab < 2 (the fab couplings 
appear with a factor of 2 in Eq. (2)) (see further remarks at the end of this section). The 
dimensional parameter fi can be constrained as follows. The Lagrangian terms in Eq. (3) 
lead to an effective quartic interaction for the h'^ and k~^~^ fields given by 



Evaluating the diagrams in Fig. 3, we obtain: 
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Since the effective couplings are negative, it follows that there must exist in the tree-level 
Lagrangian terms similar to those in Eq. (9) whose couphngs (call them A, A', A") should 
be positive and greater in absolute value than Xeff, X'^ffi Kff (otherwise the vacuum will be 
unstable). Requiring that the theory be in the perturbative regime with respect to these 
interactions (A, A', A" < 1) imposes the following constraints on the /i parameter: 



mh X (67r^)^/^ if mk ~ mh 
II < { mhX {271"^^/^ if nik < mh 
mh X (247r^)^/'^ if mk'>mh ■ 



(11) 
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FIG. 3. Diagrams giving rise to the effective quartic interaction terms for the k and h fields. 



The upper limits on the Yukawa couplings jah and hah can be made more precise by 
resorting to the renormalization group equations (RGE) and by demanding that the cou- 
plings remain perturbative to a momentum scale at least an order of magnitude higher 
than the weak scale. For example, consider the RGE evolution of the coupling f^T-: 
df^r/dt — [3/(47r^)]/^^ + ... {t = In//). Prom this it follows that an initial value of /u,- = 1 
will increase to about 1.5 at a scale an order of magnitude higher. Any larger initial value of 
f^r will signal non-perturbative effects. Similarly, the coupling Ae// of Eq. (9) will receive 
a correction proportional to /^^, obtained from the RGE dX/dt = 2/^^/7r^ + .... An initial 
value of fi^r = 1.5 will result in Ag// — 2.4 at a momentum scale an order of magnitude 
higher, which will not be in the perturbative regime. An initial value of f/^r — 1 will lead 
to Ae// = 0.47 at a scale ten times higher, which is within the perturbation theory. These 
arguments suggest that fab should be smaller than about 1 for perturbation theory to be 
valid. Similar arguments suggest that hab < 2. 



III. FITTING THE NEUTRINO OSCILLATION DATA 

In this section, we first show how the two-loop neutrino mass model fits the current 
neutrino oscillation data consistently and then derive the constraints imposed by the data 

on the model parameters. 

We use the standard parametrization of the neutrino mixing matrix (the MNS matrix) 
in terms of three angles and a phase [16]: 

/ 1 \ / ci3 sue-'^' \ / ci2 si2 \ 

C/ = C23 S23 1 -Si2 Ci2 (12) 

V -S23 C23 / V -^13e'^ Ci3 y V ^ 

where Cij = cos9ij,Sij = sinOij, and Oij is the mixing angle between the flavor eigenstates 
labeled by indices i and j. The transformation 

U^MM = Mdiag (13) 
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diagonalizes the neutrino mass matrix of Eq. (4) . 

Before going further, we shall make some remarks on the possible values of the neutrino 
masses. An interesting feature of the two-loop neutrino mass model is that, owing to the 
antisymmetric nature of the coupling matrix /, the determinant of the mass matrix M.^, is 
zero [9]. Thus one of the neutrino mass eigenvalues will be zero within the model. Although 
this result does not prevail when three-loop and higher order corrections are included, the 
induced mass at these higher loops will be extremely small and can be neglected for all 
practical purposes. 

If one uses the standard parametrization of U shown in Eq. (12), the neutrino mass eigen- 
values should be allowed to be complex in general. Since an overall phase can be removed 
by field redefinitions, there are two relative (Majorana) phases in the mass eigenvalues in 
general. However, since one of the mass eigenvalues is zero in the model under discussion, 
there is a single Majorana phase, the relative phase of the two non-zero mass eigenvalues. 
Note that with the choice of U as in Eq. (12), the elements of the matrix / cannot be made 
real. 

Experimental results indicate at least two mass hierarchies (we shall not consider the 
LSND experiment here): from the atmospheric oscillation data (z/^ ^ Ur), Am^^,„j ~ 2.5 x 
10^^ eV^ [2], and from the solar neutrinos oscillation (ug ^ or u^-) Am|j\^y^ < 10^"^ eV^ 
(the Large Mixing Angle MSW solution) or even smaller [1]. As a consequence, there are 
two possibilities for the diagonal mass matrix Aidiag in Eq. (13): the hierarchical form 

Mdiag = diag(0, m,M) with \m\ < |M| (14) 

or the inverted hierarchy form: 

M'diag± = diag(M, ±M + m, 0) with \m\ < |M| . (15) 

We have allowed two possible signs in the inverted hierarchy case. From the point of view of 
fitting the neutrino oscillation data alone, the + sign will be identical to the hierarchical form 
(Eq. (14)) since the matrices differ only by an identity matrix. However, the implications 
for the model parameters are different, and thus the predictions for rare processes will be 
different. Note also that the model does not admit the case of three-fold degenerate neutrino 
spectrum. 

The atmospheric oscillation data implies that \M\ is about 0.05 eV. This already can 
give some information on the magnitude of the parameters of the model. From Eq. (4), we 
have 

|M| « 5 X 10-12 (Ig) 



where |/| = max(/afe), \^\ = max(cja6 = mah'^i,mb). Taking / to be of order unity and /i to 
have its maximum value of about 3mh, the above equation becomes 



l/P 1^1 
rrih 



0.5 X 10-^ GeV . (17) 



Since nih has to be at least of order 100 GeV (otherwise would have already been seen 
in collider experiments), the larger of the couplings fab, uJab/GeV'^ have to be at least of the 
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order 10~^. This raises the possibility that lepton flavor violation will be observable in rare 
decays. There may even be potential conflict with the current limits. We shall discuss these 
constraints quantitatively in the next section. Our conclusion there is that presently there 
is no conflict with the data, but the rare decays should be accessible to the next round of 
experiments. 

In the rest of this section we shall perform a more detailed analysis of the neutrino 
mixings in the two-loop neutrino mass model. Our first aim is to reproduce the neutrino 
mass hierarchy and the three neutrino mixing angles. For the case of the atmospheric 
oscillations, the SuperKamiokande data indicates that the mixing is nearly maximal: |^23| — 
7r/4. On the other hand, the v^, <-> v^. mixing angle is close to zero; the reactor neutrino 
experiments [4] have set a limit sin ^13 < 0.16. For the solar neutrino oscillations, the SNO 
and SuperKamiokande results indicate that while the mixing is not exactly maximal, the 
mixing angle Q\2 is of order one, with a preferred value tan^^i2 ~ 0.4. 

In a flrst approximation let us take ^23 = 6*12 = 7r/4 and ^13 = 0. This is the exact 
bimaximal mixing limit. Then the neutrino mass matrix in the flavor basis should have the 
form 

M,^ U MdiagU'^ = - m/V2 M + m/2 M - m/2 (18) 

V M - m/2 M + m/2 / 

for the hierarchical case, or 

/ 2Af + m m /v/2 -m/v^ \ 
M',+ = U M'aiag+U^ =A M + m/2 -M - m/2 (19) 

\-m/^ -M -m/2 M + m/2 J 



Mi 



2m (-2M + m)/V2 {2M-m)/\/2' 

-2M + m)/y/2 m -m 

{2M -m)/V2 — m m 



(20) 



when the diagonal mass matrix has an inverted hierarchy form. The mass matrices above 
will get small corrections due to the deviation of the mixing angles from the exact bimaximal 
limit. 

On the other hand, by expanding Eq. (4), the neutrino mass matrix from the two loop 
model can be written in the form 



(21) 



where we have deflned 
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A 8/^ f UT J 7 * 

^ = /er/ //iT) ^ = /e/// //iT ■ (22) 

(We have written explicitly just the elements above the diagonal in Eq. (21), since the matrix 
is symmetric) . Let us assume for the moment that there is not a strong hierarchy among the 

fab, hab couplings; in this case, the hierarchy among the lepton masses {rrie ^ m^, m^) allows 
us to neglect the eUga, e'^ea quantities with respect to a;^^^, cUnr, ijJtt- With this simplification, 
the neutrino mass matrix will become 

UJrr, -^i.r (23) 

Comparing with Eq. (18), we see that a choice 

~ -^Mr ~ ^TT, e e' 1 (24) 

gives a good fit to the hierarchical neutrino mass matrix (the mass hierarchy in this case 
being m/M of order Acu/a;). For the inverse mass hierarchy case, the choice a;^^ ~ a;^^ ~ lOtt 
would fit the lower right hand side 2x2 block of the neutrino mass matrices in Eq. (19), 
but the fact that there are large elements in the first row of these matrices will require that 
there is a large hierarchy among the / couplings; that is, e, e' should be of order Mj m or 
larger. In this case, though, the simplification made to reach the form of Eq. (23) may not 
be entirely justified. 

The question arises then if it is possible to solve for the parameters e, cu exactly, without 
any approximation, in terms of the elements oi Ai^- Let's consider the general equation 

M.= ^ fuf^ = M (25) 



2 

TfJ, 



and try to solve for the parameters e, e', uJab in terms of the rriij elements of the general matrix 
Xi (the only constraints on this matrix are that it be symmetric and of zero determinant). 
At first glance, this task may seem impossible; there are five independent equations, each 
of them of order three in the unknown parameters. However, the particular form of the two 
loop mass matrix Aiu ~ f^^f'^, allows for considerable simplification of the problem. 

Note that /, having zero determinant (as a consequence of antisymmetry), has an eigen- 
vector with zero as an eigenvalue: 

<=(l,-e,e'); fvo^O. (26) 

Moreover, vq is also an eigenvector of the A4i, mass matrix; multiplying both sides of the 
matrix equation (25) by vo we obtain 

Mvo = . (27) 

The vector equation above allows solving for the parameters e, e' in terms of the elements of 
the matrix Ai (denoted by rriij): 
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mumss - mi3m23 , mi2m23 - ^1377122 
e — 2 — ' ^ — 2 — ■ \'^°) 

^22?7l33 - ^22^33 - 

The surprising result is that the e, e' parameters are completely determined by the mass 

matrix elements m.^, regardless of what the uj parameters might be. 

Once e and e' are determined, there are six 00 parameters left, and three independent 
equations. These equations are linear in uj, and therefore easy to solve. For example, we 
could use the 22, 23 and 33 matrix equahties 

C {oUrr - 2e'uJer + e'^UJee) = "^22 (29) 

C (u;^^ + 2ea;eM + ^'^^ee) = "^33 

to eliminate any three of the u parameters in terms of the other three. 

In the following, we will apply this general analysis to the two situations resulting from 
the two possible hierarchies (Eqs. (14,15)) in the neutrino masses. 

• In the case of the hierarchical neutrino mass matrix, Eqs. (28) read as 

e = tan^i2 + tan^i3 sin^23 e''^ (30) 

COS6'i3 

e' = tan^i2 ^^^^^^ — tan^i3 cos^23 e~*^ 

COS6'i3 

for general mixing angles. Taking 6'23 ~ 7r/4, ^^13 ~ 0, we have 

Then, for the LMA and LOW solutions to the solar neutrino problem, the e and e' parameters 
are in the range 0.4 — 0.5, with the difference between them being of order 2^^i3. 

As mentioned above, for the case when e, e' arc of order unity (or smaller) the euOga terms 
can be neglected with respect to o;^^, a;^^, tu,-^. Eq. (29) becomes 

+ mcos2ei2(l-^i3e'^tan^i2)/2 (32) 
C = -M/2 + m cos2^i2/2 
C ^fifM — 2 + m cos ^i2(l + ^i3e'^tan^i2)/2 

in the ^23 = 7''/4, small ^13 hmit. The relative difference between the magnitudes of the 
large uj parameters is \uj^^ — uj^T\/\^ixti\ — 2(m/M)cos^^i2, about 0.1 for the LMA solution, 
and 0.01 for the LOW solution. 

The SMA solution to the solar neutrino oscillations, although strongly disfavored by 
experimental data, can also be accommodated in this model; in this case, e, e' would be of 
order max(^i2,^i3) < 1/10. The relations (32) still hold. 

Before going further, let us consider the reverse problem, that is, determining the masses 
and mixing angles in terms of the e and uj parameters. Eq. (30) can inverted for two angles 
in terms of a third: 



10 



tan^i3 = e'%e sin^23 - e' cos^as) (33) 
tan^i2 = cos^i3(e sin^23 + 003^23) ■ 

Thus, we see that, provided that ^23 ~ 7r/4, e and e' being almost equal implies that the 
mixing angle 9is is close to zero. 

• For the inverted mass hierarchy case, the relations (28) read as 

e = -sin^23 cot^ia e~'^ , e' = cos^23 cot^ia e~'^ . (34) 

Note that here ^13 = would require f^r = 0. Barring this singular case, we see that the 

parameters e, e' are both of order O.7/613 > 5. Then, neglecting the e uj^a terms is not 
necessarily well justified. While a general analysis is possible, for the sake of simplicity we 
will restrict ourselves to the case when such terms are small. 
The equations determining the hab coupling constants are: 

M/2 + m cos ^i2(l-^i3e'^tan^i2)/2 (35) 
M/2 + m cos2^i2/2 
M/2 + m cos2^i2(l + ^i3e"'tan^i2)/2 

for the + sign in Eq. (15), and 

C corr = {m cos^^i2)/2 + M cos2^i2 (1 - 4^i3e*^tan2^i2)/2 (36) 

C ujf^r = -(m cos2^i2)/2 + (M cos2^i2)/2 

C = (m cos^^i2)/2 + M cos2^i2 (1 + 4^i3e^^tan2^i2)/2 

for the — sign. In this latter case, since m/M ~ 10^^ for the LMA solution (or smaller for the 
other solutions to the solar neutrino oscillations), and the 1-2 mixing is not quite maximal 
(cos2^i2 ~ 0.44 for LMA again), the scale in the lower right-hand corner of the neutrino 
mass matrix is set by the the terms proportional to the large mass M. This means that 
the 2-3 family symmetry visible in the exact bimaximal form of the neutrino mass matrix 
Eq. (20) docs not necessarily survive the small corrections in the mixing angles required by 
experimental results. Also, the relation (24) between the lj parameters is somewhat broken; 
we have 

Urr : u;^r : oj^^ ~ (1 - 4^i3tan2^i2) : 1 : (1 + 4^i3tan2^i2). (37) 

However, even in this case the h^^ coupling is the largest one. 

In summary, we have found that both the hierarchical neutrino spectrum and the inverted 
hierarchical spectrum with near bi-maximal mixings can be accommodated in the two-loop 
neutrino mass model without difficulty. The parameters of the model are then mostly 
determined. We now turn to other experimental signatures of the model. 

IV. EXPERIMENTAL CONSTRAINTS 

Besides neutrino masses and mixings, the Lagrangian in Eq. (1) leads to non-standard 
lepton flavor violating processes. For example, at tree level, the second line in Eq. (2) will 
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TABLE I. Constraints from lepton family number violating processes. 



allow lepton number violating decays, such as IJ~ e+e~e~ and t~ jjL~ jjL+ jjL~ , while the 
first line will give extra contributions to the standard decay of the leptons. At one loop 
level, the fab and hah couplings will contribute to the anomalous magnetic momentum of the 
e and the and will allow decays such as ^ 67. Experimental costraints will therefore 
impose limits on the parameters of the model. 

In the following, we analyze these processes in turn. 

• Lepton family number violating /i and r decays. 

The partial widths for these decay /" — > l^lc^d given by 

1 nv' 

in the limit when the masses of the decay products are neglected with respect to the mass of 
the decaying particle. The limits experimental constraints [16] set on the habhcd combinations 
are summarized in Table 1. 
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• Muonium-antimuonium oscillations. 



The last line in Table 1 is the limit arising from muonium-antimuonium conversion 
process 11+ e~ ii~e+ [17], mediated by t-channel k++ exchange. Using Fierz rearrangement, 
the amplitude for this process can be written as 

I i 

giving rise to an effective Lagrangian coupling coefficient [18] 

^MM — S 2 — ■ 

8 ml 
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• Anomalous magnetic momenta and // — > 67 -type processes. 

At the one-loop level, both /i"*" and k~^^ exchange contribute to the processes la —>■ hj- 
For the case of same leptons in the initial and final states (/„ = /(,) this leads to change in 
the anomalous magnetic moment of the lepton: 



S{9 - 2)„ 



4m| 
967r2 



iPf) 



mt 



aa {h^h)c 

mt 



For different leptons in the initial and final states, the decay width for 
rua > rrib) is: 



mt 



ml 



(39) 
Z(,7 (with 

(40) 



where a is the fine structure constant. The limits experimental results on anomalous mag- 
netic moments ^ and the /„ — > lipi processes [16] impose on the / and h couplings are 
summarized in Table 2. 

• la^ ^Jbi^b decays. 



The exchange of a h'^ boson will contribute to the semileptonic decays of the /j, and r 
leptons. Using Fierz rearrangement of spinors, the amplitude of the exchange diagram 
for the process — > I'Jb^b can be shown to be proportional to the SM diagram: 



A 



h+ 



Asm X 4 



ab 



2 



mt 



9' 



(41) 



Therefore, the angular distribution of the decay products is not affected, while the total 
decay rate is. For the case of ji decay, this imphes the redefinition of the Fermi decay 
constant (which is extracted from the muon lifetime measurements): 



2 /\/f2 



mt 



1.16639(1) X 10-^GeV" 



(42) 



where Gp refers to the SM Fermi constant, which is related to other well-measured quantities 
in the Standard Model as (in the on-shell scheme) 



Gp^ 



na 



V2M^{1-M^/Ml){l-Ar) 



(43) 



^For the anomalous magnetic moment of the muon, the bound on the non-SM contribution is 
obtained by adding the experimental error in the measurement of S{g—2)fj_ [19], the theoretical error 
in evaluating the SM prediction for this quantity [20], and the difference between the experiment 
and theory values of 6{g — 2)^ 
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Process 


Exp. bound 


Constraint 


e — > 67 


S{g - 2)e < 8. X 10-^2 


\fe^.\^ + l/erP ^ \hee\^ + M + l^erP ^ 7.6 X IQ-^ 

m| m| ^ GeV^ 




d{g - 2)^ < 1.2 X 10-8 


l/eMp + I/mtP I |/ieMp + IVl' + IVl' ^ 2.6 X IQ-^ 

m| " GeV^ 


^ 67 


Br. < 1.2 X 10-11 


\f*erUr\^ , l^ee''./- + KiJ',', + KJ> I'^r? ^ 1.7 X IQ-^^ 

ml ml ^ GeV^ 


T ^ e'-f 


Br. < 2.7 X 10-6 


l/eVMrP 1 l/l^e^er + Zi^^V + ^er^rrP ^ 2. X lO-^^ 

ml ml GeV^ 




Br. < 1.1 X 10-6 


l/eVerP ^ |/i*^/ier + ^^^^.V + ^^Ir^rrP ^8.2x10-^2 

GeV^ 






TABLE II. Constraints from — > li^j type processes. 


Here Ar encodes the effect of radiative corrections, wliich depends on the top quark and 
the Higgs boson masses. The redefinition of from Eq. (42) can be interpreted then as a 
redefinition of the Ar parameter; Ar — > Ar + 5Ar, with 



5Ar = -8^^ (44) 

Analyses indicate that variations in Ar of order 0.002 are acceptable in the Standard Model 
[21] . This would translate into a constraint on /e^ which can be found on the first line in 
Table 3. 

It is interesting to note that the exchange adds constructively to the muon decay rate. 
The neutron and the nuclear beta decay rates, on the other hand, are unchanged compared 
to the Standard Model. The value of the CKM matrix element Vud extracted from beta decay 
measurements will have to be re-interpreted in the present model. We find an upward shift 
in Vud compared to the SM by an amount given by |5Ar| of Eq. (44). Currently there is a 
2.2 sigma anomaly in the unitarity of the first row of the CKM matris: |KdP + |KsP + iKftP 
is smaller than by about 2.2 sigma [16]. With |5Ar| = 0.002, the upward shift in the value 
of Vud can nicely reconcile this anomaly within this model. 

Table 3 also contains constraints on the / couplings coming from the widths of semilep- 
tonic T decays [16] (fines 2,3). Since the redefined Fermi constant is used to compute these 
decays, the constraints are on the differences between the / couplings. Line 4 contains the 
constraint coming from e — /i universality in r decay [22] . 
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Process 


Exp. bound 


Constraint 




|<5Ar| < 0.002 


Ife^Wrnl < 1.6 X 10-8 GeV'^ 


T~ — ^ v-j-e'iJf. 


Br = 17.83 ± 0.06 % 


1 fer 1 1 fefi \ 

ml 


< 3.4 X 10-s GeV-2 




Br = 17.37 ± 0.07 % 


|//ir| ~ |/e/i| 


< 4. X 10-8 GeV-2 


e/fj, universality 


= 0.999 ± 0.003 


|/er| ~ Iffirl 

ml 


< 2.5 X 10-8 GeV-2 



TABLE III. Constraints from /j, decay and semileptonic r decays. 



• Neutrinoless double beta decay. 



From the approximate forms of the neutrino mass matrix given in Eqs. (18)- (20), it 
follows that the effective neutrino mass relevant for neutrinoless double beta decay is ap- 
proximately m(/3/30z/) ~ m/2 ~ 10^"^ eV in the hierarchical well as in the inverted 
hierarchical case of Eq. (20). This will be difficult to observe in the near future. On the other 
hand, in the inverted mass hierarchy of Eq. (19), the effective mass is m{l3p0u) ~ M ~ 0.05 
eV, which should be observable in the next round of experiments. 

Now, let us look at these constraints in light of the relations among the fab and hab 
couplings imposed by neutrino masses and mixings results of Sec. III. Consider first the / 
parameters. In the case of hierarchical form for the neutrino mass matrix, from Eqs. (30), 
(31) we have 

fefj, ~ fer ~ ■ (45) 

The strongest constraint on the / couplings comes then from the // — > 67 process. Prom 
Table 2 we have 



Jj. 
2m? 



< 0.4 X 10-8 GeV-2 . (46) 



On the other hand, from the last of Eq. (32), we have 



(167r2)2 ml 2 
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with M ~ -\/Z5~xT(P^ eV = 5 X 10 ^^GeV (here we have neglected the contribution of the 
terms proportional to mg) . Plugging in the numbers we find 

rrih ~ 10^ (-^) fl^K^ i GeV . (48) 

Assume for now that is smaUer or of the same order as mh', then I k, 2. Using the upper 
limits on the h couphngs and n/nih ratio derived in Sect. II, we get the following allowable 
range for the /i+ scalar mass: 

10"^ f^r GeV < nih < lO' GeV . (49) 

Here the lower limit comes from the // — > 67 constraint Eq. (46), while the upper limit 
comes from the neutrino mass equations above. We see then that when the /^^ couphng 
takes its maximum value (of order unity), the h'^ scalar will be out of reach of future colliders. 
However, this is the most unfavorable case; if the h, /, and couplings take values smaller 
than the highest values admissible, the upper bound on mh moves to lower values. For 
example, if /^r ~ 0.1, ruh has to be smaller than 10 TeV. Note also that Eqs. (46), (48) 
impose a lower limit on the strength of the /^,-, /i^^ couplings: 

hfxiji > 

This requires that ruh has to be greater than about 100 GeV. Also, the lower bound on h^^ 
justifies neglecting the Wge, f-^e/* terms in the neutrino mass matrix; indeed, for the ujen term 
to be of the same order of magnitude as uj^n, we would need he^ ~ 200 x h^^ which will 
move the theory into the nonperturbative regime. 

Wc conclude from the preceding discussions that the decay ^ — e+7 should be accessible 
to the next round of rare decay experiments. Some of the model parameters are already 
excluded by current limits. For example, if nih ~ 1 TeV, any value of /^^ > 2 x 10~^ will be 
inconsistent with current limits. There are plans to improve the present limit on — > e + 7 
by several orders of magnitude in the near future [24] . This decay is predicted to be within 
reach of these improved experiments. 

The leptonic phenomenology constraints on the hab coupling constants and the mass of 
the doubly charged scalar ruk are somewhat weaker. All constraints from processes involving 
the electron can be made to go away by choosing hee, and/or her to be close to zero (which 
is allowed by the neutrino masses and mixings pattern). Then, the strongest constraint on 
the hab/rrik ratio comes from the t~ — > process; using Eq. (24) to relate hf^r and 

h^^, this constraint reads as 

^ < 1.3 X 10-=^ GeV-^ 
rrik 

which can be easily satisfied without violating the other constraints discussed previously. 
Note also that it is not required that hee, h^fx and her be close to zero; independent constrains 



0.1 
0.1 



far I \^hj 



— J >l-xlO-^ 
mhj 

-^^ ' > 1.7 X 10-2 



(50) 
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FIG. 4. Scalar boson production cross-section 



on these quantities from Table 1 and Table 2 are: 



\h h* I 
ml 



< 



4. X 10 



-7 



ml 



< 



4. X 10"^ 



GeV^ 



ml 



< 



9. X 10"^ 



GeV^ 



allowing any of them to be of the same order of magnitude as /i^^ (although not both hee 
and hen the same time). 

The bounds on the parameters of this model are more stringent for the inverted mass 
hierarchy case. In this case |/e^| ~ |/er| — |/^r|/sin^i3 are the larger / couphngs. In terms 
of fefx the bounds on m^ can be written: 



2Vsin^i3 X 10^ /e^ GeV < mn - 2sm^9n x 10^ f^^h^^ I GeV 



(51) 



< 2sin2^i3 X 10^ fl GeV 



with fen < 1. The first inequality comes from the ^ 67 process in Table II, while the last 
one comes from the neutrino mass equations (35) (we considered the + case here). Note that 
the allowed interval is narrower and shifted toward lower values of mh- The lower bounds on 
the fe/x and and /i^^^ in Eqs. (50) are increased by a factor l/(sin^i3)^/^. This also imposes 
a lower bound on the value of the ^13 mixing angle: 



sin^i3 > 0.046 



(52) 



which should be testable in neutrino oscillation experiments. These constraints are relaxed 
by a factor cos2^i2(l + 4^i3tan2^i2) (~ 0.44(1 + 8^13) for tan^^is = 0.4) for the - sign in 
the inverted mass hierarchy case (using Eqs. (36) for the neutrino mass matrix elements). 
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V. COLLIDER SIGNALS 



We now turn to the possible signals coming from the production of the charged scalars 
at colliders. Consider first the case of hadron colliders. Since the h'^ and k'^'^ scalars do 
not couple directly to the quarks, the production of these particles will proceed through the 
s-channel processes: qq j*,Z* h~^h~,k'^'^k . The cross-section for these processes, 
at the 2 TeV Tevatron and at the LHC, are presented in Fig. 4, as a function of the mass 
of these particles. As can be seen from Fig. 4, the cross-section for the production of the 
doubly charged scalar k'^'^ is about four times larger than the one for the production of the 
singly charged scalar /i+ (for equal masses). Experimentally, the is also much easier 
to see; the hierarchy between the h couplings hfj,fj, : h^r '■ — 1 : Trifx/'mr '■ {fnn/rrirY 
(Eqs. (24,37)) imphcs that the will decay predominantly to a same sign muon pair (or 
electron pair, if h^e is of the same order of magnitude as h^j^fj). In any case, the experimental 
signature of 4-lepton final state will be striking. The SM background from ZZ production 
can be greatly reduced by imposing high pt cuts on the transverse momenta of the leptons, 
and by requiring that the invariant mass of opposite sign lepton pairs be different from the 
Z mass. In Ref. [23], it has been estimated that as few as 10 events are enough for the 
discovery of a doubly charged boson which decays mostly to e or pairs. In this case, the 
Tevatron Run II will be able to probe up to about 250 GeV with an integrated luminosity 
of 15 fb-\ while the LHC reach will be about 800 GeV with 100 fb-\ or 1 TeV with 1 ab-^ 
integrated luminosity. 

At a hadron collider, the s-channel production of a h^h~ pair will be much harder to 
detect experimentally. The final state of two leptons -|- missing energy (associated with 
the neutrinos coming from decay) has SM backgrounds coming from ZZ, ZW (with one 
lepton lost) and WW production, as well as Drell-Yan production of two leptons with mis- 
measured energy, cuts can be used to reduce the background in this case too, but a 
more detailed analysis is needed to obtain the hadron collider reach for production in 
this mode. Another interesting possibility would be the production of a single through 
radiation from the v line in the process qq' — > lui. The final state in this case will be three 
leptons plus missing Et- High pt cuts and cuts on the invariant mass of lepton pairs can be 
used here too to reduce the SM background (which will come mostly from WZ production). 
However, the production cross section in this case is proportional to the / couplings squared. 
For boson masses below 1 TeV, the ji ej constraint Eq. (46) requires that / be in the 
0.01 - 0.1 range; as a consequence, there will be at most of order tens of single events 
produced at the LHC for any /i+ mass in the TeV range. 

At e~^e~ colliders, s-channel pair production of singly charged or doubly charged bosons 
will be limited by the energy of the machine. However, due to the cleanliness of the envi- 
ronment, they will be easy to sec, if kincmatically accessible. The reach in mass at these 
machines will therefore be roughly y/s/2. Provided that the h^e coupling is large, k~^~^ can 
also be singly produced at an e~e~ collider. More interesting, though, would be the study 
of the doubly charged scalar at a muon collider, since the /i^^ coupling of the k'^'^ has to be 
large in this model. 
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VI. CONCLUSIONS 



In this paper we have performed a detailed analysis of a specific two-loop neutrino 
mass model. We have shown that the model can accommodate both the solar and the 
atmospheric neutrino data, while at the same time satisfying all bounds arising from leptonic 
phenomenology. While at present there is no conflict with any hmit, the model predicts that 
the rare decays r ^ 3fi and /i — e + 7 should be within reach of forthcoming experiments. 

Although the model contains two independent coupling matrices, one symmetric and 
the other antisymmetric, they are well constrained by the current neutrino oscillation data. 
We have found a simple way of relating these matrices to neutrino oscillation parameters. 
In particular, we have shown that the relative ratios between the three elements of the 
antisymmetric coupling matrix / are completely determined by the neutrino mixing angles, 
regardless of the other parameters of the model. Since the contribution of the symmetric 
couphngs hab to the neutrino mass matrix is proportional to the product of the masses of 
the a and b charged leptons, only /i^^, /i^^ and /i,-,- are relevant for neutrino oscillations. The 
facts that the atmospheric mixing angle 623 is close to 7r/4 and that there is a hierarchy 
in the solar and atmospheric neutrino mass-splittings fixes the relative ratio between these 
three parameters to be : h^r '■ hrr — 1 : ^n/rrir : {mn/mry. The fact that /i^^ is the 
larger of these may have implications on the search for the doubly charged scalar A;"*""*" in 
collider experiments. 

From the experimental result on the atmospheric neutrino mass splitting and the con- 
straints on the anomalous process fi 67, a lower bound of 0.01 on the largest of the / 
and h coupling constants has been derived. There are upper bounds on the masses of the 
charged scalars h'^ and k~^~^ as well. For non-maximal values of the Yukawa couplings (for 
example, f^r 0.1) the predicted values for the scalar masses are in the TeV range, 

which should be probed at the LHC, and perhaps at Run II of the Tevatron. 
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